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Abstract: We follow the evolution of an asymptotically AdS black brane with a fixed 
temperature gradient at spatial infinity until a steady state is formed. The resulting energy 
density and energy flux of the steady state in the boundary theory are compared to a con¬ 
jecture on the behavior of steady states in conformal field theories. Very good agreement is 
found. 
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1 Introduction 


Predicting the out of equilibrium behavior of many-body systems is a notoriously difficult 
problem. Many of the tools which have been developed to describe thermodynamic equilib¬ 
rium are insufficient to explore systems which are far from it. Of the various non-equilibrium 
phenomenon, one that seems to have a tractable lead is the steady-state problem where one 
attempts to study the properties of a steady state which emerges once the system is placed 
between two heat baths. Recent progress on this problem has been seen in a variety of sys¬ 
tems which may be described by a quantum field theory, see, e.g., [1—8]. In this work we focus 
on the steady-state generated when placing a 3 space-time dimensional relativistic conformal 
gauge-theory between two heat baths. 

When placing a system described by a two dimensional conformal field theory between 
two heat baths its late time, steady state, behavior is completely fixed in terms of the left and 
right central charges, the levels of the chiral currents (if present) and the relative pressure 
difference of the heat baths [1-5]. A recent conjecture which appeared in [4-6] attempts to 
generalize this result to higher dimensional theories. 

The form of the conjecture of [4-6] is as follows. Let us consider two heat baths which are 
infinitely separated along one spatial direction, call it z, and are translation invariant along 
the remaining directions, x±. Let us denote the pressure of the left and right heat baths by 
Pl and Pr respectively, and the relative pressure difference by 


dp 


Pl-Pr 
Pl + Pr 


( 1 . 1 ) 


Without loss of generality we choose Pl > Pr. Consider an initial state whose energy 
momentum tensor smoothly interpolates between the two heat baths and is also translation 
invariant along x±. The mechanism by which energy is propagated in such a setting can 
be categorized into two types. A flow driven (or ballistic) mechanism by which the energy 
propagates via sound modes or their non-linear counterpart, and diffusion. If the system is 
flow driven then we may expect that at late times the energy of the system will propagate 
with a velocity vl towards the left heat bath and with velocity vr towards the right heat 
bath generating a steady-state in between. If, in addition, the system thermalizes sufficiently 
rapidly so that the disturbance immediately behind the propagating waves is in thermal 
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equilibrium, then energy-momentum conservation fixes the expectation value of the pressure, 
T zz , and the energy flux along the 2 direction, T tz , at late times [4-6]. 

More precisely, in [5] it was shown that the aforementioned physical considerations about 
the evolution of the system leads to two branches of solutions for the possible late time 
behavior of T zz and T tz . In [4, 6] it was argued that a translation invariant steady state 
allows for only one possible late time solution. The work of [5] and [4, 6] are consistent. 

In 3 space-time dimensions, the translation invariant steady-state (referred to as the 
thermodynamic branch in [5]) is characterized by 


= ^ (4 - \fl-5p- 


P , TV (pT z 

—— = lim -—— 

P() t—>00 Pq 

Y = Urn Tr p T ' * = 

Tq t->oo Jt’q 3 v 


( 1 . 2 ) 


where Pq is the average pressure 


^0 = ^ (Pl + Pr) , 


(1.3) 


and p represents the density matrix for the steady state. The expressions for J and P in the 
other branch of the steady state solution is rather long and has been omitted. A plot of J 
and P as a function of 5p for both branches appears in figure 1. The reader is referred to [5] 
for details. 
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Figure 1. A plot of the conjectured universal values of the late time steady state energy flux (left 
panel) and pressure (right panel). Here P 0 is the average pressure of the heat baths, given by (1.3), and 
Sp is the relative pressure difference, given by (1.1). The red curves correspond to the thermodynamic 
branch of steady-states where the expectation value of the energy momentum tensor is equivalent to 
that of a boosted thermally equilibrated configuration. The blue curve corresponds to the “other” 
branch of steady-state configurations whose physical role is somewhat obscure. 


The late time steady-state prediction described in (1.2) has been verified explicitly for 
systems whose time evolution is determined by hydrodynamics [4-6] . While the hydrodynamic 











description is probably a good one at small relative pressure differences, dp <C 1, it is unlikely 
to be valid when the relative pressure difference of the heat baths is large and there are large 
pressure gradients when one moves away from the steady state towards null infinity. In this 
work we compare the conjecture of [4-6] to steady states in three dimensional theories which 
may be described using the AdS/CFT correspondence, up to 5p = 0.5, and find excellent 
agreement. 


2 Constructing black brane steady states 


Simply stated, the AdS/CFT correspondence [9-11], or more generally gauge-gravity duality 
(see, e.g., [12]), relates a classical d + 1-dimensional gravitational theory to a planar, d- 
dimensional conformal held theory (CFT) at strong coupling. Of particular importance to 
the problem at hand is that solutions to classical gravity which follow from the action 


s = ^I^( R + it) dix (21) 

and asymptote to conformally hat space at spatial infinity hold information about a dual 
conformal held theory in 3 space-time dimensions. 

For instance, the black brane solution to the equations of motion which follow from (2.1) 
after setting LAdS = 1 is given by 


ds 2 


2 dtdr — r 




( 2 . 2 ) 


with T the Hawking temperature of the black brane. The solution (2.2) is dual to a thermal 
state given by a density matrix p = of the conformal held theory, with H its Hamiltonian 
and /3 = T^ 1 its inverse temperature. The expectation value of the energy momentum tensor 
of the CFT in this thermal state is given by 


TV ( pT ***') 


/ 2P(T) 0 0 \ 

0 P(T) 0 

^ 0 0 P(T) J 


(2.3) 


with 

P(T) = Po( 4 f^y (2.4) 

where the value of the dimensionless parameter po depends on the particular details of the 
dual theory under consideration. In the (planar and strongly coupled) ABJM theory [13], 


2N 2 

P0 ~972X 


(2.5) 


with N the rank of the gauge group and A = N/k with k the Chern-Simons level. 
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In order to generate a steady state in the dual CFT of the type to which the conjecture 
of [4-6] applies, we consider black branes whose dual energy momentum tensor will satisfy 


and 


(2P(T r ) 0 0 \ 

jiin Tr(pT flu (t,x ± ,z)) = 0 P(T R ) 0 

infinite \ 0 0 -P(Tr) J 


lim Tr (pT^(t,x±, z)) 

z—>■—oo 

t,Xj_ finite 


/ 2P(T l ) 0 0 \ 

0 P(T L ) 0 
v 0 o P(T L )J 


(2.6a) 


(2.6b) 


with T r and Tj j the temperatures of the right and left heat baths respectively and p the 
density matrix of the system. Geometrically, the conditions (2.6) correspond to a metric 
which asymptotes to (2.2) as z approaches ±oo with corresponding temperature Tr and Tr. 
An initial configuration for the metric which interpolates between the two asymptotic black 
brane solutions will evolve in time and a time independent configuration is expected to be 
generated at late times and finite 2 . Once such a steady state is made available then the 
energy momentum tensor associated with the late time geometry of such a black brane can 
be compared to the prediction of [4-6]. This is the strategy we will follow in the remainder 
of this work. Other examples of out of equilibrium stationary black hole configurations can 
be found in [14-17]. 

Consider the metric ansatz 


ds 2 = 2 dt (dr - A(t , r)dt - F(t, z, r)dz ) + S 2 (t, *, r) (e B(t ’ z ' r) dx\ + e~ B ^ z ’ r) dz 2 ^j . (2.7) 

As pointed out in [18] this ansatz is invariant under r —> r + £(t, z). We will fix this residual 
diffeomorphism symmetry shortly. The equations of motion which follow from (2.1) take the 
form 


4d 2 £ + £ ( d r B f = 0 (2.8a) 

-T 2 d 2 r F - T 2 d r B d r F + C F [B, £]T = 5 F [S, B] (2.8b) 

4£ 3 3 r £ + 4£ 2 a r ££ = S’ s [£, B, F] (2.8c) 

A^drB + 4£ 3 5 r £B = 5^[£, B , F, ±] (2.8d) 

2£ 4 5 2 A = 5a[S, B , F, E, B] (2.8e) 

and 

E = Q^[E, B, F, A] (2.9a) 

d z F = Q dzF [T, B, F, A] , (2.9b) 

with the following definitions. Dotted variables are given by 

X = d t X + Ad r X X = d t X + Ad r X . (2.10) 
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The expressions Cf and Sx in (2.8) depend only on spatial derivatives of their arguments, 
and the expressions Qx in (2.9) depend on spatial and time derivative of their arguments. 
The asymptotic AdS boundary, located at r —> oo, is given by the line element 


ds 2 = ( 2 dtdr + r 2 (— dt 2 + dx]_ + dz 2 )) (l + O (r 4 )) . 


( 2 . 11 ) 


The energy momentum tensor of the dual theory can be read off of the large r asymptotics 
of (2.7) as follows. Given (2.11) we find that, under the equations of motion, the large r 
expansion of the metric components A, F, T, and B are given by 

A = \ (r + £(t, z)) 2 - dt$(t, z) + ° i r ~ 2 ) 


_ „ . fi(t,z) 3d z bAt,z) „ , ox 

F = -d z £(t, z) + \ , + 77 ; A 9 + 0 i r ) 

r + £(t,z) 4 (r + £(t,z)) 2 v ' 

S = r + z) ~ 40(r 3 + ( ^^ ))5 + ° K 6 ) 


( 2 . 12 ) 


B = 


h (t,z) 


+ O (r” 4 ) . 


(r + C(t,z)) 3 

Here a\(t,z), f\(t,z) and b^(t,z) must satisfy the constraints 

d t ai = |9 2 /i d t fi = ^d z ai + d z b 3 


(2.13) 


(which are a result of (2.9)), but are otherwise undetermined. Once (2.13) are satisfied then 
(2.8) ensure that (2.9) will hold. Following [19, 20] the energy momentum tensor of the dual 
theory is given by 


/—2ai |/i 


Tr(pTH=Po 


0 

f/i -ai - §63 0 

0 0 -ai + ^b 3j 


(2.14) 


where p is the density matrix associated with the dynamical black brane solution. In order 
to generate the steady state ( 2 . 6 ) we will impose 


p 0 ai ( t , z 00) = - P(T r ) Poai ( t , z -s> -00) = - P(T L ) 

and also the initial conditions 

/i(t = 0 , z — >• ±00) = 0 b 3 (t = 0 , z — > ±00) = 0 . 


(2.15) 


(2.16) 


which ensure, by causality, that f\ and 63 vanish at large \z\ for any t. 

Our strategy for solving (2.8) is identical to that presented in [18]. At t = 0 we set 


B = 0 


fi =0 


a\ = —Aq ( 1 — a tanh (/3 tanh ( — 


(2.17) 
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Using (2.17) and fixing the residual diffeomorphism symmetry parameter £(t) we can now 
determine E, F, E, B and A at t = 0 by solving (2.8) in sequence together with the asymptotic 
conditions (2.12). From this data, one can determine B , a i and f\ at t = At using (2.10) 
and (2.13). This procedure can now be repeated to obtain the fields at t = 2At and so on, 
ad infinitum. 

The numerical scheme which was implemented in order to evolve the metric forward in 
time was based on pseudo-spectral methods to solve the linear equations (2.8) and 3rd order 
Adams-Bashforth to evolve forward in time. We used a grid of 23 points along the radial 
direction, r, and 100 points along the z direction which we parametrically compactified using 

C = tanh(^/(AL)) (2-18) 

for an optimal value of L. Rescaling all the boundary coordinates by a factor of A is a 
symmetry of the equations of motion due to conformal invariance of the boundary theory. In 
practice we have implemented a code which uses A = 1. In solving the equations of motion 
we found it convenient to subtract from A, B , E, E and F the solution to (2.8) at t = 0 with 
the initial conditions (2.17) which can be obtained analytically, viz., 

S = r + C F = -d g £ A= l -{r + tf-d t { i + ^- 

2 r + £ 

S = ^(r + 0 2 + ^ B = 0. (2.19) 

In practice the value of £ was chosen so that an apparent horizon will be located at r = 1 
which is the endpoint of our grid. We refer the reader to [18] for details of such a scheme. 

The final expression for the energy momentum tensor of the dual CFT has been computed 
using (2.14). We have carried out the computation for several values of 0 < 5p < 0.7. A 
typical result is plotted in figure 2. For 5p < 0.5 we have managed to evolve our numerics 
for sufficiently long so that the energy flux approached its equilibrium value to an accuracy 
of less than 1 %. See figure 3. At very large relative pressure difference the black brane 
configuration at z —> oo has a very small temperature relative to the one on at z —> — oo and 
our numerical scheme seems to break down. At best this breakdown is due to the parametric 
compactihcation we have used to map the z coordinate to a finite interval. In order to 
avoid difficulties associated with such compactihcations [ 21 ], it may be possible to repose the 
steady state problem using periodic boundary conditions in the z direction. However, if the 
numerical instability associated with large relative pressure differences is due to the formation 
of caustics, then one would need a more robust algorithm than that of [18] and summarized 
here to obtain a handle on the problem. 

When dp 2 ~ 0.9 the steady state pressure and energy density of the other (blue) branch 
becomes lower than that of the thermodynamic (red) branch. It would be interesting to check 
if at such large pressure differences the resulting steady state prefers the other branch. 
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Figure 2. A typical plot of the energy flux (left) and pressure (right) of a theory dual to the black 
hole described by (2.7). These plots were generated by solving (2.8) with initial conditions given in 
(2.17) with A = 1, A 0 = 100, a ~ 0.8657 and /3 = 1/2 (corresponding to 5p = 0.4). At late times the 
steady state value of the energy flux and pressure asymptote to the red branch to within 0.2%. 

Max(J/P 0 ) AJ 



Figure 3. Plots describing the approach to the conjectured steady state configuration. These plots 
were generated by solving (2.8) with initial conditions given in (2.17) with A = 1, Aq = 100, a ~ 1.082 
and /? = 1/2 (corresponding to 5p = 0.5). The z coordinate was parametrically compactified as in 
(2.18) with L = 5. In the left panel we have plotted, for every time t, the maximal (in z) value for 
the energy flux J together with the conjectured steady state value of J (red and blue). In the right 
panel we have plotted the relative difference between the maximal (in z) value for the energy flux J 
and the predicted values of the thermodynamic (red) and other (blue) branch. 
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